library(readxl)
Belgium <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Belgium.xlsx")
View(Belgium)
# Importing Data
BELGIUM <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Belgium.xlsx")
# Creating Time Series Data
BELGIUM_ts <- ts(BELGIUM, start=c(2004,1), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
BELGIUM_ts
sum(is.na(BELGIUM_ts))
library(forecast)
BELGIUM_ts <- tsclean(BELGIUM_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(BELGIUM_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
BELGIUM_ts_model <- auto.arima(BELGIUM_ts)
BELGIUM_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
BELGIUM_ts_forecast <- forecast (BELGIUM_ts_model, level=c(95), h=258)
plot(BELGIUM_ts_forecast)
BELGIUM_ts_forecast
write.table(BELGIUM_ts_forecast, file="Belgium_TSA.csv", sep=",")
# Importing Data
CANADA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Canada.xlsx")
# Creating Time Series Data
CANADA_ts <- ts(CANADA, start=c(2007,1), end=c(2019,03), frequency=12)
# Viewing and Checking the Created Time Series Data
CANADA_ts
sum(is.na(CANADA_ts))
library(forecast)
CANADA_ts <- tsclean(CANADA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(CANADA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
CANADA_ts_model <- auto.arima(CANADA_ts)
CANADA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
CANADA_ts_forecast <- forecast (CANADA_ts_model, level=c(95), h=261)
plot(CANADA_ts_forecast)
CANADA_ts_forecast
write.table(CANADA_ts_forecast, file="Canada_TSA.csv", sep=",")
# Importing Data
AUSTRALIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Australia.xlsx")
# Creating Time Series Data
AUSTRALIA_ts <- ts(AUSTRALIA, start=c(2004,1), end=c(2017,10), frequency=12)
# Viewing and Checking the Created Time Series Data
AUSTRALIA_ts
sum(is.na(AUSTRALIA_ts))
library(forecast)
AUSTRALIA_ts <- tsclean(AUSTRALIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(AUSTRALIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
AUSTRALIA_ts_model <- auto.arima(AUSTRALIA_ts)
AUSTRALIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
AUSTRALIA_ts_forecast <- forecast (AUSTRALIA_ts_model, level=c(95), h=278)
plot(AUSTRALIA_ts_forecast)
AUSTRALIA_ts_forecast
write.table(AUSTRALIA_ts_forecast, file="Australia_TSA.csv", sep=",")
# Importing Data
INDONESIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Indonesia.xlsx")
# Creating Time Series Data
INDONESIA_ts <- ts(INDONESIA, start=c(2006,01), end=c(2019,01), frequency=12)
# Viewing and Checking the Created Time Series Data
INDONESIA_ts
sum(is.na(INDONESIA_ts))
library(forecast)
INDONESIA_ts <- tsclean(INDONESIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(INDONESIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
INDONESIA_ts_model <- auto.arima(INDONESIA_ts)
INDONESIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
INDONESIA_ts_forecast <- forecast (INDONESIA_ts_model, level=c(95), h=263)
plot(INDONESIA_ts_forecast)
INDONESIA_ts_forecast
write.table(INDONESIA_ts_forecast, file="Indonesia_TSA.csv", sep=",")
# Importing Data
FRANCE <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/France.xlsx")
# Creating Time Series Data
FRANCE_ts <- ts(FRANCE, start=c(2004,1), end=c(2019,01), frequency=12)
# Viewing and Checking the Created Time Series Data
FRANCE_ts
sum(is.na(FRANCE_ts))
library(forecast)
FRANCE_ts <- tsclean(FRANCE_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(FRANCE_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
FRANCE_ts_model <- auto.arima(FRANCE_ts)
FRANCE_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
FRANCE_ts_forecast <- forecast (FRANCE_ts_model, level=c(95), h=263)
plot(FRANCE_ts_forecast)
FRANCE_ts_forecast
write.table(FRANCE_ts_forecast, file="France_TSA.csv", sep=",")
# Importing Data
INDIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/India.xlsx")
# Creating Time Series Data
INDIA_ts <- ts(INDIA, start=c(2007,04), end=c(2019,04), frequency=12)
# Viewing and Checking the Created Time Series Data
INDIA_ts
sum(is.na(INDIA_ts))
library(forecast)
INDIA_ts <- tsclean(INDIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(INDIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
INDIA_ts_model <- auto.arima(INDIA_ts)
INDIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
INDIA_ts_forecast <- forecast (INDIA_ts_model, level=c(95), h=260)
plot(INDIA_ts_forecast)
INDIA_ts_forecast
write.table(INDIA_ts_forecast, file="India_TSA.csv", sep=",")
# Importing Data
GERMANY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Germany.xlsx")
# Creating Time Series Data
GERMANY_ts <- ts(GERMANY, start=c(2004,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
GERMANY_ts
sum(is.na(GERMANY_ts))
library(forecast)
GERMANY_ts <- tsclean(GERMANY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(GERMANY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
GERMANY_ts_model <- auto.arima(GERMANY_ts)
GERMANY_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
GERMANY_ts_forecast <- forecast (GERMANY_ts_model, level=c(95), h=257)
plot(GERMANY_ts_forecast)
GERMANY_ts_forecast
write.table(GERMANY_ts_forecast, file="Germany_TSA.csv", sep=",")
# Importing Data
CHINA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/China.xlsx")
# Creating Time Series Data
CHINA_ts <- ts(CHINA, start=c(2005,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
CHINA_ts
sum(is.na(CHINA_ts))
library(forecast)
CHINA_ts <- tsclean(CHINA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(CHINA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
CHINA_ts_model <- auto.arima(CHINA_ts)
CHINA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
CHINA_ts_forecast <- forecast (CHINA_ts_model, level=c(95), h=257)
plot(CHINA_ts_forecast)
CHINA_ts_forecast
write.table(CHINA_ts_forecast, file="China_TSA.csv", sep=",")
# Importing Data
AUSTRIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Austria.xlsx")
# Creating Time Series Data
AUSTRIA_ts <- ts(AUSTRIA, start=c(2004,1), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
AUSTRIA_ts
sum(is.na(AUSTRIA_ts))
library(forecast)
AUSTRIA_ts <- tsclean(AUSTRIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(AUSTRIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
AUSTRIA_ts_model <- auto.arima(AUSTRIA_ts)
AUSTRIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
AUSTRIA_ts_forecast <- forecast (AUSTRIA_ts_model, level=c(95), h=258)
plot(AUSTRIA_ts_forecast)
AUSTRIA_ts_forecast
write.table(AUSTRIA_ts_forecast, file="Austria_TSA.csv", sep=",")
# Importing Data
USA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/USA.xlsx")
# Creating Time Series Data
USA_ts <- ts(USA, start=c(2005,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
USA_ts
sum(is.na(USA_ts))
library(forecast)
USA_ts <- tsclean(USA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(USA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
USA_ts_model <- auto.arima(USA_ts)
USA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
USA_ts_forecast <- forecast (USA_ts_model, level=c(95), h=257)
plot(USA_ts_forecast)
USA_ts_forecast
write.table(USA_ts_forecast, file="USA_TSA.csv", sep=",")
# Importing Data
KUWAIT <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/KUWAIT.xlsx")
# Creating Time Series Data
KUWAIT_ts <- ts(KUWAIT, start=c(2008,01), end=c(2018,12), frequency=12)
# Viewing and Checking the Created Time Series Data
KUWAIT_ts
sum(is.na(KUWAIT_ts))
library(forecast)
KUWAIT_ts <- tsclean(KUWAIT_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(KUWAIT_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
KUWAIT_ts_model <- auto.arima(KUWAIT_ts)
KUWAIT_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
KUWAIT_ts_forecast <- forecast (KUWAIT_ts_model, level=c(95), h=264)
plot(KUWAIT_ts_forecast)
KUWAIT_ts_forecast
write.table(KUWAIT_ts_forecast, file="Kuwait_TSA.csv", sep=",")
# Importing Data
UK <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/UK.xlsx")
# Creating Time Series Data
UK_ts <- ts(UK, start=c(2004,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
UK_ts
sum(is.na(UK_ts))
library(forecast)
UK_ts <- tsclean(UK_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(UK_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
UK_ts_model <- auto.arima(UK_ts)
UK_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
UK_ts_forecast <- forecast (UK_ts_model, level=c(95), h=257)
plot(UK_ts_forecast)
UK_ts_forecast
write.table(UK_ts_forecast, file="UK_TSA.csv", sep=",")
# Importing Data
TAIWAN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Taiwan.xlsx")
# Creating Time Series Data
TAIWAN_ts <- ts(TAIWAN, start=c(2002,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
TAIWAN_ts
sum(is.na(TAIWAN_ts))
library(forecast)
TAIWAN_ts <- tsclean(TAIWAN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(TAIWAN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
TAIWAN_ts_model <- auto.arima(TAIWAN_ts)
TAIWAN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
TAIWAN_ts_forecast <- forecast (TAIWAN_ts_model, level=c(95), h=258)
plot(TAIWAN_ts_forecast)
TAIWAN_ts_forecast
write.table(TAIWAN_ts_forecast, file="Taiwan_TSA.csv", sep=",")
# Importing Data
ITALY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Italy.xlsx")
# Creating Time Series Data
ITALY_ts <- ts(ITALY, start=c(2004,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
ITALY_ts
sum(is.na(ITALY_ts))
library(forecast)
ITALY_ts <- tsclean(ITALY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(ITALY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
ITALY_ts_model <- auto.arima(ITALY_ts)
ITALY_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
ITALY_ts_forecast <- forecast (ITALY_ts_model, level=c(95), h=258)
plot(ITALY_ts_forecast)
ITALY_ts_forecast
write.table(ITALY_ts_forecast, file="Italy_TSA.csv", sep=",")
# Importing Data
SPAIN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Spain.xlsx")
# Creating Time Series Data
SPAIN_ts <- ts(SPAIN, start=c(2004,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
SPAIN_ts
sum(is.na(SPAIN_ts))
library(forecast)
SPAIN_ts <- tsclean(SPAIN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(SPAIN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
SPAIN_ts_model <- auto.arima(SPAIN_ts)
SPAIN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
SPAIN_ts_forecast <- forecast (SPAIN_ts_model, level=c(95), h=258)
plot(SPAIN_ts_forecast)
SPAIN_ts_forecast
write.table(SPAIN_ts_forecast, file="Spain_TSA.csv", sep=",")
# Importing Data
THAILAND <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Thailand.xlsx")
# Creating Time Series Data
THAILAND_ts <- ts(THAILAND, start=c(2007,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
THAILAND_ts
sum(is.na(THAILAND_ts))
library(forecast)
THAILAND_ts <- tsclean(THAILAND_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(THAILAND_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
THAILAND_ts_model <- auto.arima(THAILAND_ts)
THAILAND_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
THAILAND_ts_forecast <- forecast (THAILAND_ts_model, level=c(95), h=257)
plot(THAILAND_ts_forecast)
THAILAND_ts_forecast
write.table(THAILAND_ts_forecast, file="Thailand_TSA.csv", sep=",")
# Importing Data
KOREA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Korea.xlsx")
# Creating Time Series Data
KOREA_ts <- ts(KOREA, start=c(2008,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
KOREA_ts
sum(is.na(KOREA_ts))
library(forecast)
KOREA_ts <- tsclean(KOREA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(KOREA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
KOREA_ts_model <- auto.arima(KOREA_ts)
KOREA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
KOREA_ts_forecast <- forecast (KOREA_ts_model, level=c(95), h=257)
plot(KOREA_ts_forecast)
KOREA_ts_forecast
write.table(KOREA_ts_forecast, file="Korea_TSA.csv", sep=",")
# Importing Data
JAPAN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Japan.xlsx")
# Creating Time Series Data
JAPAN_ts <- ts(JAPAN, start=c(2001,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
JAPAN_ts
sum(is.na(JAPAN_ts))
library(forecast)
JAPAN_ts <- tsclean(JAPAN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(JAPAN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
JAPAN_ts_model <- auto.arima(JAPAN_ts)
JAPAN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
JAPAN_ts_forecast <- forecast (JAPAN_ts_model, level=c(95), h=258)
plot(JAPAN_ts_forecast)
JAPAN_ts_forecast
write.table(JAPAN_ts_forecast, file="Japan_TSA.csv", sep=",")
# Importing Data
SLOVAKIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Slovakia.xlsx")
# Creating Time Series Data
SLOVAKIA_ts <- ts(SLOVAKIA, start=c(2004,05), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
SLOVAKIA_ts
sum(is.na(SLOVAKIA_ts))
library(forecast)
SLOVAKIA_ts <- tsclean(SLOVAKIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(SLOVAKIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
SLOVAKIA_ts_model <- auto.arima(SLOVAKIA_ts)
SLOVAKIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
SLOVAKIA_ts_forecast <- forecast (SLOVAKIA_ts_model, level=c(95), h=258)
plot(SLOVAKIA_ts_forecast)
SLOVAKIA_ts_forecast
write.table(SLOVAKIA_ts_forecast, file="Slovakia_TSA.csv", sep=",")
library(readxl)
Taiwan <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Taiwan.xlsx")
View(Taiwan)
# Importing Data
TAIWAN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Taiwan.xlsx")
# Creating Time Series Data
TAIWAN_ts <- ts(TAIWAN, start=c(2002,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
TAIWAN_ts
sum(is.na(TAIWAN_ts))
library(forecast)
TAIWAN_ts <- tsclean(TAIWAN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(TAIWAN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
TAIWAN_ts_model <- auto.arima(TAIWAN_ts)
TAIWAN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
TAIWAN_ts_forecast <- forecast (TAIWAN_ts_model, level=c(95), h=258)
plot(TAIWAN_ts_forecast)
TAIWAN_ts_forecast
write.table(TAIWAN_ts_forecast, file="Taiwan_TSA.csv", sep=",")
# Importing Data
TAIWAN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Saudi Arabia/Taiwan.xlsx")
# Creating Time Series Data
TAIWAN_ts <- ts(TAIWAN, start=c(2002,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
TAIWAN_ts
sum(is.na(TAIWAN_ts))
library(forecast)
TAIWAN_ts <- tsclean(TAIWAN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(TAIWAN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
TAIWAN_ts_model <- auto.arima(TAIWAN_ts)
TAIWAN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
TAIWAN_ts_forecast <- forecast (TAIWAN_ts_model, level=c(95), h=258)
plot(TAIWAN_ts_forecast)
TAIWAN_ts_forecast
write.table(TAIWAN_ts_forecast, file="Taiwan_TSA.csv", sep=",")
